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ABSTRACT 



Context. The first unidentified very high energy gamma ray source (TeV J2032+4130) in the Cygnus region has been the subject of 
intensive search for a counterpart source at other wavelengths. In particular, observations in radio and X-rays are important to trace a 
population of non-thermal electrons. 

Aims. A deep (~ 50 ksec) exposure of TeV J2032+4130 with XMM-Newton has been obtained. The large collection area and the 
field of view of the X-ray telescopes on-board of XMM-Newton allow to search for faint extended X-ray emission possibly linked to 
TeV J2032+4130. 

Methods. The contribution of point sources to the observed X-ray emission from TeV J2032+4130 is subtracted from the data. The 
point-source subtracted X-ray data are analyzed using blank sky exposures and regions adjacent to the position of TeV J2032+4130 
in the field of view covered by the XMM-Newton telescopes to search for diffuse X-ray emission. 

Results. An extended X-ray emission region with a full width half maximum (FWHM) size of as 12 arc min is found. The centroid 
of the emission is co-located with the position of TeV J2032+4130. The angular extension of the X-ray emission region is slightly 
smaller than the angular size of TeV J2032+4130 (FWHM=14 ± 3 arc min). The energy spectrum of the emission coinciding with the 
position and extension of TeV J2032+4130 can be modeled by a power-law model with a photon index T = 1.5 ±0.2 stat ±0.3 sys and an 
energy flux integrated between 2 and 10 keV of /2-10 kcv ~ 7 • 10~ 13 ergs/(cm 2 s) which is lower than the very high energy gamma-ray 
flux observed from TeV J2032+4130. The energy flux detected from the extended emission region is about a factor of two smaller 
than the summed contribution of the point sources present. The energy spectrum can also be fit with a thermal emission model from 
an ionized plasma with a temperature k B T « 10 keV. 

Conclusions. We conclude that the faint extended X-ray emission discovered in this observation is the X-ray counterpart of 
TeV J2032+4130. Formally, it can not be excluded that the extended emission is due to an unrelated population of faint, hot 
(k B T ss 10 keV) unresolved point-sources which by chance coincides with the position and extension of TeV J2032+4130. We 
discuss our findings in the frame of both hadronic and leptonic gamma-ray production scenarios. 



1. Introduction 

Ground based air Cherenkov telescopes have detected in the last 
6 years a number of unidentified sources of very high energy 
(VHE, E > 100 GeV) gamma-rays located in the Galactic plane. 
The first representative of these sources (TeV J2032+4130) was 
discovered serendipitou sly with the HEGRA telescope system 
in the Cygnus region dRowell & Hornsl 120021 : lAharonian et al.L 
120021) . 

In a subsequent deep exposure for a total of * 280 hours 
of observation time the source detect ion was confirm e d and 
its properties were stu died in detail dHorns & Rowe 1 I2004t 
lAharonian et al.Ll2005al) . 

A remarkable feature of TeV J2032+4130 is the fact that 
the angular extension of the source is (6.2 + 1.2 stat + 0.9 sys ) 
arc min in radius which corresponds to a FWHM=14 arc min. 
Gamma-ray emission from TeV J2032+4130 has been con- 
firmed independently by the W hipple collaboration dLang et all 
2004; Konope lko et all l2006l) . A gamma-ray excess from the 



* Based on observations obtained with XMM-Newton, an ESA sci- 
ence mission with instruments and contributions directly funded by 
ESA Member States and NASA. 



Cygnus region albeit mor e extended was re cently reported by 
the Milagro collaboration dAbdo et all 120061) . The properties of 
TeV J2032+4130 are common to most of the other unidenti- 
fied VHE gamma-ray sources that have been discovered with 
the H.E.S.S. telescopes so far dAharonian et all l2005bl l2006al) : 
TeV J2032+4130 can be considered as the proto-type of the 
unidentified VHE gamma-ray sources known today. 
So far, a number of point-like or moderately extended candi- 
dates for cou nterparts of TeV 12032+4130 have b een identi- 
fied at radio dParedes et all l2006tlButt et alll2006bl) an d X-rays 
dButt et alll2003ll2006atlMuk heriee et al., 2003i|2006|). 
The origin and nature of TeV 12032+4130 remains unclear. It 
has been suggested that the nearby (d = 1.7 kpc) massiv e 
stellar OB association Cyg OB2 (see e.g. iKnodlsedCT (|2000)) 
could be an accelerator of charged c osmic rays and conse- 
quently a s ite of gamma-ray produc t ion (lAharonian et all |2002; 
iButt et all l2006at lAharonia n et all l2006bl) . This scenario has 
been strengthened by the recently discovered spatially extended 
VHE gamma-ray emi ssion from the directio n of the open stellar 
cluster Westerlund 2 dAharonian et allfeOOTl) . 
Other possible explanations for the nature of TeV 12032+4130 
have been brought forward including gamma-ray production in 
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possible jet lobes of Cyg X-3 ( Aharonian et al.L 120021) . an un- 
known pulsar wind nebula (Bednarekl 120031), or even extra- 
galact ic s ource candidates a s suggested by iMukheriee et alj 
d2003h and lButt et alJ (l2006bl) . 

Independent of the nature of the source, two different gamma- 
ray production mechanisms are generally considered. In the 
hadronic scenario, gamma-rays are produced mainly via in- 
elastic scatterin g of accelerated nuclei with the ambien t 
medium (see e.g. | Aharonian et al.l f2005al) : lReimer et all d2006): 
iDomingo-Santamaria & Torres! (|2006)). Alternatively, in the 
leptonic scenario gamma-rays are produced via inverse Compton 
scattering of ambi ent photons from a population of ener - 
getic electrons (lAharonian et all l2005at iReimer et all 2006). 
Recently, it has been suggested, that excitation of giant dipole 
resonances of relativistic heavy nuclei in radiation dominated 
environme nts could be responsible fo r gamma-ray production in 
Cyg OB2 dAnchordoqui et alll2006l) . 

While the gamma-ray observations so far have been inconclusive 
with respect to the origin of the observed signal, X-ray observa- 
tions can provide additional information to identify the origin of 
the emission and to discern between the two proposed scenarios. 
In this Letter, we report the detection of spatially extended X-ray 
emission coinciding with the position of TeV J2032+4130. 

2. Observations and Data analyses 

The data were taken during two separate pointings of XMM- 
Newton. Table Q] summarizes the available data and configura- 
tions of the instruments used. The data were screened for soft 
proton flares following the method suggested in lRead & Po nman 
(2003). The detectors performed almost nominally: CCD #6 of 
the MOS 1 camera had been switched off and CCD #5 of MOS2 
shows an increased instrumental background below 1 keV which 
is two times larger than the background seen in the adjacent 
CCDs #4 and #6. However, above 1 keV the background rate in 
CCD #5 appeared consistent with the other detectors. The data 
suffer from contamination of single scattering events from the 
bright X-ray source Cyg X-3. The contamination is most promi- 
nent in the energy band above 5 keV. For this reason, the analy- 
ses that have been performed are constrained on the energy range 
from 1-5 keV where the instrumental background and stray light 
contamination from Cyg X-3 are minimal. A closer inspection of 
the contamination using ray-tracing simulations has shown that 
the contribution of scattered light from Cyg X-3 to the X-ray 
emission seen from TeV J2032+4130 is negligible. At the same 
time, t he vignetting of the telescopes is sm allest in this energy 
range (lKirschLl2006tlPradas & Kerpll2005l) . 
Data reduction and analyses were performed using the Standard 
XMM-Newton Science Analysis Software SAS v7.0. 

2.1. Search for diffuse X-ray emission 

In an initial step of the data analyses, a catalogue of sources 
detected in the 1-5 keV energy range is assembled using the 
edetect_chain task simultaneously for the MOS1, MOS2, and 
EPIC pn cameras for each of the two exposures. The minimum 
detection likelihood is set to a low value of 2. This way, 13.5 % 
of the detected sources are due to statistical fluctuations while at 
the same time, more and fainter sources are detected than with 
the default value of 10. A detailed study of the point sources de- 
tected in this XMM-Newton observation will be presented else- 
where. Here, the source catalogue has been used to generate 
model images of the point sources separately for each camera 
and pointing taking into account the position dependent point 



spread function of the XMM-Newton mirrors. These model im- 
ages are then subtracted off the individual frames of the three 
cameras and for the two pointings. Within the 6.2' region cov- 
ered by the extension of the TeV source, a total energy flux of 
/i-5 keV ~ 10~ 12 ergs/(cm 2 s) associated with point sources is 
observed. In order to subtract the particle and extragalactic pho- 
ton background blank sky event files for the corresponding ob- 
servation modes and filter settings have been obtained from the 
XMM guest observers' facility^]. The blank field data have been 
processed to match the observational data and the resulting im- 
ages are subtracted off the pointings towards TeV J2032+4130. 
The images for the MOS 1 and MOS2 cameras obtained from the 
two exposures are combined to form a background subtracted 
mosaic. The images of the EPIC pn camera have been omit- 
ted because of artefacts resulting from over-subtracting the con- 
tribution of point sources. For each of the four remaining im- 
ages, exposure maps have been generated using the SAS task 
eexpmap. The exposure maps are combined in a mosaic which 
is then trimmed to a minimum exposure of 5 ksec corresponding 
to roughly 5 % of the peak value. The background subtracted 
mosaic is then divided by the exposure map to obtain an expo- 
sure corrected, background subtracted, and source free image. 
Fig. Q] shows three images of the combined exposures of the 
MOS 1 and MOS2 cameras before the subtraction of the model 
images of the point sources (a), after the subtraction and smooth- 
ing with Gaussian of fixed width of 45 arc sec (b), and finally, 
the image after convolving it with a Gaussian of fixed width of 
1.5 arc min (c). Note the presence of scattered light along ring 
segments in the southern part of Fig. la. While Fig. lb shows 
indications for the presence of a r im-like feature resem bling the 
morphology observed in radio bv lParedes et alj d2006l) . it should 
be noted that the subtraction of point-sources has an impact on 
the observed morphology on angular scales small compared to 
the extension of the detected extended emission. For a compari- 
son with the extension and morphology of the gamma-ray source 
TeV J2032+4130, contours are overlaid to Fig. QJ. The X-ray 
image shows an extended emission region which is co-located 
with TeV J2032+4130 and of similar extension. The exposure 
corrected image of the diffuse emission is fit by a Gaussian 
with a width of (5.1 + 0.3) arc min which corresponds to a 
FWHM= (11.7 ± 0.6) arc mirfl 



2.2. Energy spectrum of the diffuse X-ray emission 

In order to determine the energy spectrum of the diffuse X-ray 
emission, the contribution of the point sources has to be removed 
by excising regions around the sources. We use an excising re- 
gion adapted to the point spread function and brightness of the 
source. The size of the region is chosen such that the relative 
contribution of the source to the local background is less than 
20 % at the boundary of the excised region. The wings of the 
point spread function lead to a contamination of the detected dif- 
fuse emission which amounts to 27 % of the excess signal seen. 
This effect of contamination would be however compensated by 
the opposite effect caused by the diffuse emission present in the 
excised region (» 25 % ) and therefore excluded from the re- 
constructed energy spectrum. We estimate the systematic uncer- 
tainty of this contamination to be 10 % of the total flux. The en- 
ergy spectra of a source region with a radius of 6.2' centered on 
the position of TeV J2032+4130 from the MOS2 and EPIC pn 



1 ftp://xmiti.esac.esa.int 

2 all errors given throughout the paper are to the confidence of one 
standard deviation 
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Fig. 1. Background subtracted and exposure corrected images obtained from the combined exposures of the MOS1 and MOS2 
cameras (in units of counts/(sec arcmin 2 ). While the left image shows the emission including all point sources, the middle image 
shows the image after subtracting off all point sources and smoothed with a Gaussion of 45 arc sec width. Finally, the right image 
is obtained by smoothing the source subtracted image with a Gaussian of 1 .5 arc min width and changing the scaling to highlight 
the extended X-ray emission present. The green contours indicate (in linear spacing) the significance contour of the HEGRA 
observations (starting at 3 cr). 



cameras of the second pointing (ObsID 0305560201) are ex- 
tracted. The second pointing is centered 5' north of the centroid 
of TeV J2032+4130 and therefore, we can choose a background 
region which is mirrored through the center of the field of view 
and has a similar acceptance as the source extraction region. 
Using the background region in the same field of view, the back- 
ground estimate includes all relevant background components 
including X-ray emission from the Galactic ridge and the par- 
ticle and instrumental background specific to this observation. 
The small azimuthal modulation of the acceptance of the MOS 
cameras can be neglected along this direction. The MOS1 cam- 
era is not used because the disabled CCD#6 unfortunately does 
not allow to cover the entire source. 

The response files are generated for these particular observations 
using the standard tools arfgen and rmfgen. Finally, the ge- 
ometrical areas of the source and background regions are cal- 
culated taking gaps between the CCDs, bad pixels, and point 
source exclusion regions into account. The resulting background 
subtracted X-ray spectrum is fit with a power-law as well as an 
optically thin hot plasma model (apec) including photoelectric 
absorption (phabs). The spectral fitting of the two spe ctra was 
done using the xspec vl 1.3.2p spectral fitting package ( Arnaud, 
1996). The fit describes the data well (see Table[2]for a summary 
of the fit parameters; only statistical 1 <x errors are quoted). 
In order to estimate the influence of faint, unresolved point 
sources present in the stellar cluster, we have reduced X-ray 
imaging data taken with the Chandra satellite (Obs.# 4501 for 
48.6 ksec on TeV J2032+4130 and Obs.# 4511 for 98.7 ksec 
in the south-west of TeV J2032+4130, centered on the core of 
Cyg OB2). Based upon a comparison of the Chandra and XMM- 
Newton source catalogues and the study of the X-ray source pop - 
ulation present in the Cyg OB2 cluster (Colom bo et al.L 12006). 
we estimate that point sources below the XMM-Newton detec- 
tion limit can contribute up to « 30 % of the observed diffuse 
emission. It is however not straight-forward to correct the ob- 
served excess for the unresolved point sources detected with 
Chandra as the flux of the sources is generally found to vary 
with time. We estimate the systematic uncertainty on the flux 
conservatively to be 50 % combining the influence of the unre- 



Table 1. Summary of observation times and configurations on 
the target "TeV J2032+4130". All instruments were operated in 
full frame mode with a medium filter. The exposure quoted in 
parenthesis is for the EPIC pn camera. 



ObsID 


Date 


R.A. 


Dec. 


Exposure 








J2000 


J2000 


[ksec] 


0305560101 


2005- 


10-21 


20:31:57 


41:29:58 


27.3(23.6) 


0305560201 


2005- 


10-25 


20:31:57 


41:34:55 


25.8(20.5) 



Table 2. Summary of the fit results of the X-ray energy spectrum 
with a power-law model (powerlaw) or thermal emission (apec 
with fixed solar abundances of the plasma) including photo- 
electric absorption (phabs). 



Parameter 


Value (powerlaw) 


Value (apec) 


N H [10- 1 cm--] 


3.5 ± 1.6 


3.2+1.1 


T, k B T [keV] 


1.5 ±0.2 


10.5 + 3.2 


/2-10 keV [10~ 13 ergs/(cm 2 s)] 


7.3 ± 1.1 


7.4+1.1 


XL (d-o-f.) 


1.02 (862) 


1.01 (861) 



solved point sources (a 30 %), uncertainties on the background 
subtraction (« 10 %), as well as the tails of the point spread func- 
tions of excised point sources (as 10 %). The total, unabsorbed 
energy flux derived from the model fit (see Table |2]i is found to 



be/ : 



2-10 



(7.3 + 3.5 sys + l.l stat ) x 10-° ergs/(cm 2 s) 



3. Conclusions 

The XMM-Newton observations presented here indicate the 
presence of an extended (FWHM=11.7 arc min) X-ray source 
co-located with the first unidentified VHE gamma-ray source 
TeV J2032+4130 discovered with the HEGRA air Cherenkov 
telescope system. The size of the X-ray source is similar to the 
one of TeV J2032+4130. The energy spectrum can be fit by a 
power-law model or by a thermal emission model with a plasma 
temperature of kgT * 10 keV. The unabsorbed energy flux of 
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E[eV] 

Fig. 2. The spectral energy distribution of TeV J2032+4130: the 
bow-tie indicates the energy spectrum as measured with XMM- 
Newton including systematic and statistical uncertainties. The 
lines indicate the result of a model calculation for a hadronic 
gamma-ray and neutrino production scenario for a continuously 
active accelerator of protons up to a maximum energy of 5 PeV 
at an age of f age = 2 500 yrs and magnetic field B = 1 mG. 
The radio flux point of a possible non-therm al extended radio- 
source at A = 20 cm dParedes et al.L 120061) is assumed to be 
an upper limit to the actual radio emission associated with 
TeV J2032+4130. The upper limits between 0.5-5 keV (CX, 
Chandra) and 20- 40 keV (ISGRI, INTEGRAL) are taken from 
iButt et ail d2006al) . The energy spectrum of the EGRET source 
3EG J2030+4118 dHartmanet all 1 19991) is considered to be an 
upper limit. The Milagro (MI) upper limit i s the integrated emis- 
sion in a 3x3 (°) 2 region ( Abd o et alll20 06) while the lower limit 
(MI) is scaled to the solid angle covered by TeV J2032+4130 
assuming a uniform surface brightness. The upper limit from 
the HEGRA- AIROBICC (HA) wide angle Cherenkov detector 
is taken from lAharonian et~ai1 d2002l) . 



- — — Young source (1500 yrs) ' 3egj2033+4hs' 




E[eV] 

Fig. 3. In the leptonic scenario, the data are described well by a 
young source of an age of 1500 yrs, a magnetic field of B — 3 fiG 
and an energy density of the seed photon field of wir = 3 eV/cm 3 
with a grey body temperature of T = 10 K. As an alternative, an 
older source at an age of 35 kyrs is shown (B = 3 fiG, Wm 
I eV/cm 3 ). For a description of the multi-wavelength data see 
the caption of Fig. [2] 



the X-ray source in the energy range from 2-10 keV is a factor 
of 2-3 smaller than the one observed from TeV J2032+4130 at 
energies from 1-10 TeV. 

We note that the observed extended X-ray emission could in 
principle be a so far unknown population of faint X-ray sources 
that are by chance distributed at an angular size which is sim- 
ilar to the one of TeV J2032+4130. The energy spectrum of 
these sources would have to be different from the bulk of the 
stella r X-ray sources detected from Cyg OB2 (IColombo et all 
2006) which on average show a thermal spectrum with kgT ~ 
1 ... 3 keV. Taking these considerations into account, we con- 
clude that the observed extended X-ray emission is the X-ray 
counterpart of TeV J2032+4130. 

Initial modeling of the X-ray and gamma-ray energy spectra us- 
ing a hadronic gamma-ray production scenario indicates that the 
observed spectra can be naturally explained by a young (a few 
kyrs) "Pevatron" accelerator. The observation of synchrotron X- 
ray emission up to » 5 keV constrains the product of the square 
of maximum energy of accelerated protons (E^) and magnetic 
field (B) to exceed £^ ax • B ^ 5 PeV 2 mG. In this picture, the 
gamma-ray energy spectra would continue without cut-off well 
beyond 10 TeV. The hard X-ray emission from a similar source 
size region as the gamma-ray emission is a natural prediction 
within this model where the X-rays are produced by synchrotron 
emission of secondary electrons. Notably, in the frame of this 
model, the extended gamma-ray emission found by the Milagro 
collaboration would be explained by gamma-rays produced by 
accelerated particles of energies ^100 TeV which have already 
escaped the accelerator. 

A leptonic scenario provides a valid explanation of the observa- 
tions as well (see Fig. |3J- In this scenario, a young accelerator 
with a low magnetic field of a few /iG accelerating electrons fol- 
lowing a power-law distribution with dN/dE oc E~ 2 reaching up 
to energies of a few 100 TeV can provide a good fit to the X- 
ray and gamma-ray data. In contrary to the hadronic scenario, 
the gamma-ray energy spectrum above 10 TeV is expected to 
be rather soft due to un-avoidable Klein-Nishina suppression of 
inverse Compton scattering. Furthermore, the seed photon den- 
sity has to be higher than the average value in the interstellar 
medium (which is not unlikely given the possible proximity to 
the Cyg OB2 stellar association and the large stellar extinction 
towards that region indicating high density of dust). In the lep- 
tonic scenario, the X-ray spectrum can be expected to vary at 
different parts of the source due to cooling effects. More X-ray 
observations will be required to detect spectral variability. 
Finally, neutrino observations will prove decisive to discern be- 
tween the two emission scenarios (since neutrinos are only ex- 
pected in the hadronic scenario). Even though the neutrino flux 
calculated here is very likely not detectable wit h the coming gen- 
eratio n of neutrino telescopes like IceCube dAchterberg et aU 
2006), it may be that the gamma-ray flux detected with Milagro 
from the same region is tracing the high energy particles which 
have already left the accelerator and fill a much larger volume. In 
this case, the region around TeV J2032+4130 would be a pow- 
erful high energy neutrino source and would be dete ctable with 
future neutrino telescopes (Beacom & Kistler[ |2007l) . 
It should be noted that the observations of TeV J2032+4130 car- 
ried out with XMM-Newton are so far the deepest observations 
available for any of the unidentified VHE gamma-ray sources. 
It will be interesting to perform similar observations to search 
for faint, extended X-ray counterparts of other unidentified VHE 
gamma-ray sources. 
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